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^ , Abstract. 

Mh' Strong evidence for pairing and superfluidity has recently been found in atomic 

^ ' Fermi gases at the BCS-BEC crossover both in cohective modes and RF excitation 

energies. It is argued that the scale for the effective pairing gaps measured in RF 
experiments is set by the lowest quasiparticle in-gap excitation energies. These are 
calculated at the BCS-BEC crossover from semiclassical solutions to the Bogoliubov- 
i-rt ! deCennes equations. The strong damping of the radial breathing mode observed in 

^ ' the BCS limit occur when the lowest quasiparticle excitation energies coincide with 

the radial frequency, which indicates that a coupling between them take place. 
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tJ- ■ 1. Introduction 

t^ , 

\Q , The recent developments in cold Fermi gases mark another milestone a century 

2 ' of marvelous discoveries within the field of superfluidity and superconductivity. 

O '• Experiments have established the stability of cold Fermi gases with strong attractive 

'^ ! interactions PJ El El IH Ej and proven the unitarity limit near Feshbach resonances. A 

smooth BCS-BEC crossover is found in which the Fermi atoms gradually bind into 
^ , (Bose) molecules near Feshbach resonances as predicted in Refs. [El El IE] ■ The collective 

C ! modes provide strong evidence for superfluidity jHl 1101 ■ RF spectroscopy ^^ clearly 

Q ! show a peak in the response function that depends on interaction strength, density and 

temperature as expected from resonance superfluidity theory [T^IT!^. 
k> ! The advantage of the trapped atomic gases is that we can tune the interactions, 

densities, temperatures and other trap parameters in a controlled way and thus 
explore superfluidity in detail as well as in general. This insight has already been 
exploited to describe some features of pairing in nuclei and neutron stars [HI ITKj . 
Furthermore, the pairing gaps A are of order the Fermi energy Ep around the BCS- 
BEC crossover and therefore A/Ep is an order of magnitude larger than in high 
temperature superconductors and superfluid ^He, and two orders of magnitude larger 
than in standard superconductivity (see, e.g, [1^1) ■ Just a few years ago such large 
pairing gaps was not considered possible by most of the condensed matter community. 
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In a uniform system the gap is A = O.bAEp in the unitarity hmit according to Monte 
Carlo calculations J^. The observable transition frequencies are, however, related to 
the quasiparticle energies which are considerably smaller than the pairing field in the 
center of the trap and therefore in-gap excitations. 

The purpose of this work is to compare the recent experimental results on RF 
spectroscopy and collective modes with theoretical calculations. After a general 
introduction to pairing and superfluidity in uniform systems in section II, we address 
pairing in traps in section III with special emphasis on the in-gap quasiparticle 
excitations at the BCS-BEC crossover, and how they compare to the data of ^J in 
section IV. In section V we discuss the collective mode frequencies and damping and 
the evidence for superfluidity, and finally give a summary. 

2. Pairing in Uniform systems 

The quasiparticle excitations in traps are usually very different from those in uniform 
system. Yet the bulk pairing field is important for understanding the BCS-BEC 
crossover, and it is very useful as a basis for pairing in traps with many particles within 
the Thomas- Fermi approximation (TF). 

Solving the gap equation at zero temperature for a Fermi gas interacting through 
an attractive scattering length a < gives a pairing gap in the dilute limit, akp -C — 1, 

IT 

_2akF 

Here, k = 8/e^ in standard BCS. Gorkov included induced interactions which reduces 
the gap by a factor ~ 2.2 to k = {2/ey^^ fTj. In the unitarity limit kpla]^!, which 
can be reached around the Feshbach resonances, the gap is of order the Fermi energy 
[3 El E]. Extrapolating (0) to akp -^ ±oo gives a number close to that found from 
odd-even staggering binding energies A = OAdEp calculated by Monte Carlo |14j . 

The pairing gap can qualitatively be followed in the crossover model of Leggett [Zj 
by solving the gap equation 

i = ^!^y[l-±l. (2) 



A = kEp exp 



(1) 



Here, the quasiparticle energy is E-^ = \/(£k — yu)^ + A^ with e^ = h'^k'^/2m. The 
chemical potential fi follows from conservation of particle number density 

^ = Efi-^l- (3) 

k L Ek } 

In the dilute (BCS) limit the gap equation leads to the standard BCS gap of Eq. ((H) - 
not including the Gorkov correction (see Fig. 1). The chemical potential is /i = Ep and 
does not include the standard mean field Hartree-Fock correction 27ran/m of a dilute 
gas. 

The advantage of the crossover model is that it extends to the strongly interacting 
(molecular BEG) limit. Here, the pairing gap approaches A = AEp/y/Snakp. The 



chemical potential approaches half of the molecular binding energy E^ = —Ti^ /ma? with 
a BEC mean field corresponding to a molecular scattering length of a^ = 2a. Four- 
body jIHI, Monte Carlo calculations [TH] and experiments jlj do, however, indicate that 
am — 0.6a 

On the BCS side [akp < 0) the minimum quasiparticle energy is A and occur when 
k = kp- On the BEC side [akp > 0) the chemical potential is negative and the minimum 
quasiparticle excitation energy is the quasiparticle energy for k = 



E = Jfi^ + A^. (4) 



The chemical potential is observed in the spin excitation response function 

Monte Carlo calculations [HI Hn] of binding energies, equation of states and pairing 
gaps agree qualitatively with the crossover model. The pairing gap is A/Ep = 0.54 in 
the unitarity limit [H] which is quite close to the extrapolation of the Gorkov gap 
^Gorkov I j^^ _ {2/ey^'^ ~ 0.49 but somewhat smaller than that of the crossover model 

^Leggettj^^ = 0.69. 

In the BCS limit the critical temperature for the BCS transition is T^ = 
(e*"/7r)A, where C = 0.577.. is Eulers constant. In crossover model there are 
two transition temperatures in the BEC limit jH], namely the molecular BEC one 

at Tc = [n/2^(3/2)]^/^7r/?7i = 0.218-E'i;' and the molecule dissociation temperature 
rpdissoc ^ Eb/HEb/Epf/^ 

3. Pairing in harmonic oscillator traps 

The various regions of pairing were described in pl| for a dilute system of A^ atoms 
with Hamiltonian 

^ a 

H=Y: Hoir,) + Ann'- ^ ^^(r, - r,) , (5) 

i=i "^ i<j 

in a harmonic oscillator (h.o.) potential Ho{r) = p^/2m + J2l=i'i^'-^k''^k/'^- We shall 
mention a few relevant results only before we investigate the strongly interacting limit. 

At least two dimensionless parameters are required to describe this system even in 
the spherically symmetric case: ui = 002 = uj^ = ujq. These are, e.g., the number of 
particles N and the interaction strength a (when energies are measured in units of huo 
and lengths in aosc = d'h/mwo). Several interesting pairing regions or "phases" appear 
vs. N and a. In contrast, the binding energies and pairing in an uniform gas at zero 
temperature are functions of one parameter only, e.g. akp. 

When the traps contain relatively few atoms, iV^SlO'^, that are weakly interacting, 
the mean field does not significantly lift the degenerate angular momentum states 
/ = np^ Up — 2, ..., 1 or 0, where np = Ep/fiuj = {3NY^^, due to the SU(3) symmetry of 
the spherical symmetric h.o. potential. Consequently, pairing takes place between all 
these states which leads to the supergap ^Q 



32y/2np \a\ 32 

A = G ^ ^^^—^-0 = j^M-inu^o . (6) 



Here, kp = \/2np/aosc — l-^A^^^^/ctosc is the Fermi wave number in the center of the 
trap. For more particles in the trap the stronger mean field cause level splitting, which 
reduce pairing toward single level pairing, A^^^^, which displays a distinct shell structure 
with h.o. magic numbers. Pairing in nuclei has a similar shell structure on top of an 
average gap equal to the supergap which for constant density scales with the atomic 
mass number A as A ^ G ~ 5.5MeV /A^'^. For very large nuclei pairing approaches 
that in bulk matter A ~ l.lMeV. 

For stronger interactions pairing also takes place between shells and the gap 
increases to [21] 

" l-21n(e^n^)G/^cuo' ^ ^ 

This multi-shell pairing region extends up to 2G\n{e'"nF)^huj. 

For even stronger interactions the pairing field exceeds the harmonic oscillator 
energy in the center of the trap and the coherence length, C, = /cf/ttttiA, is much 
smaller the TF radius Rtf of the cloud. This spatially inhomogeneous case with a 
strong pairing field A(r) can be solved by the Bogoliubov-deGennes equations 

E^M^(r) = [Ho + f/(r) - fi]ur,{r) + A(r)t;^(r) , 
Er,v^{r) = -[Ho + f/(r) - /i]t;^(r) + A{r)u^{r) . 

Here, E^, are the quasiparticle energies, u^ and f^ the Bogoliubov wave functions, and 
U{r) the mean field. The pairing field can be approximated by the Gorkov gap in the 
TF approximation |21j . 

Let us first study the spherical symmetric trap. Tthe Bogoliubov wave function can 
be written on the form m^ = Uni{r)Yim{9,(j)), where (/,m) are the angular momentum 
quantum numbers. The Bogoliubov-deGennes equations can be solved semiclassically 
|22j which leads to a WKB quantization condition 



(n+ -)- = m dr , ^ (8) 

2 2 Jr, ^2m(/i - f/(r)) - ryat,, - l{l + l)/r^ 

Here, n is the number of radial nodes which in the dilute limit is the h.o. number 
counted from the Fermi level. R12 are the classical turning points. When the pairing 
field is strong, i.e. much larger than the quasiparticle excitation energies, it determines 
the inner turning point by A(i?i) = £"„;. Ri is then close to the outer turning point 
i?2 — Rtf- The single particle excitations therefore take place near the surface of the 
trap where the pairing field is weak. Furthermore, since R2 — Rtf the centrifugal 
potential can be neglected when the angular momentum is small, I -^ np, and we shall 
denote this set of energies by En- 

In the dilute BCS limit /i ~ Ep and R2 = Rtf = V^^Fdosc in Eq- dHJ, and the 
quantization condition reduces to 



1 7C [Rtf jEl-A{rY 



i-TF 



independent of / as long as / -C nj?. In this expression we can use TF for the 



pairing field A{kp{r)) of Eq. (jT)), i.e. with kplr) = J2nF{i — r"^ /Rtf)/(^osc and 
Ef{t) = h kp(r)/2m. The resulting quasiparticle energies En i21j are given by inverting 
the relation 



kF\a\huo ~ - — " < In 
(n + 1/2) [ 



K,{n + l/2)^7r^ni7 {huo) 



-1 



(10) 



4 E^ 

n 

It is valid above the multishell pairing regime of Eq. (|7j) and up to the dense limit 
l/{akF)< - 1. 

We can furthermore extend the semiclassical model to stronger interactions and the 
unitarity limit. As in the Leggett crossover model we assume that the correlations are 
included in the pairing field and that mean field can be included in /i and in an effective 
h.o. potential. The Bogoliubov-deGennes equations can then be applied - specifically the 
semiclassical quantization condition of Eq. (jH)). However, although the crossover model 
provides a simple and qualitatively description of the ground state around the BCS-BCS 
crossover, it differs from the ground state found by Monte Carlo calculations and 4-body 
calculations for both the pairing field and the chemical potential as mentioned above. 
We will therefore include the effect of mean field by replacing the chemical potential by 
that calculated by Monte Carlo (TUdnj. Induced interactions we include by using the 
Gorkov result of Eq. (^ for the pairing field which also is a fair approximation around 
the unitarity limit according to Monte-Carlo calculations |14j . 

Eq. (jHI) can now be solved in the Thomas-Fermi approximation. For example, 
in the unitarity limit A(r) = K,EF{r) and the attractive mean field is U = (3EF{r), 
where j3 = Eint/E^in is an universal parameter in the unitarity limit [13] and is directly 
related to the chemical potential at zero temperature. It has been measured directly 
from expansion energies Ql 12111 El and confirms Monte Carlo calculations (3 ~ —0.56 in 
the unitarity limit [T^ ITU] . The integral in eq. (jHl) can be expanded for radii near the 
surface and calculated analytically 



E,^ 



Lr(7/4)^ 



'^' (^'^^^)'^' n 11 



Note that the form of the mean field in the unitarity limit is such that it can be 
included in the chemical potential and h.o. potential and "renormalizes" these quantities 
by factors of {1 + f3) to some power. The mean field decreases the TF trap size 
Rtf oc {1 + py^^, and the chemical potential but increases the central density. As result 
the quasiparticle energies increase slightly. In Eq. (fTTj) the Fermi energy is defined as 
Ef = uf^uq without mean field corrections as also defined in the experiment of [TT] . 

Also note that the finite size of the system is manifest because E^/Ef is not only 
a function of akF but also depends on uf = EF/hu. For example. En/ Ef scales with 
particle number as ~ N'"^'^ in the unitarity limit. 

Around the unitarity limit the semiclassical quantization condition is solved 
numerically with the TF Gorkov gap and a mean field correction (3 averaged over 
densities in the trap. On the EEC side the chemical potential [i should be included in 
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Figure 1. Pairing gaps in units of the Fermi energy vs. interaction strength akp 
in the center of the trap. The full and dashed lines are the lowest two quasiparticle 
excitation energies Eq and Ei (see text). The "effective" pairing gaps hi'/2 measured 
in PJ at three Fermi energies are shown with symbols. 



excitations as in Eq. (j3)). For the lowest quasiparticle excitations near the surface region, 
where the density is low, this chemical potential is within TF simply the molecular 
binding energy, /i = Eb/2 = —h /2'ma?. The quasiparticle energy En is therefore 
replaced by 



' E'^ + (i?b/2)2. The lowest quasiparticle energies Eq and Ei including this 
molecular binding energy correction are shown in Fig. 1 as function of l/{akp). They 
reduce to Eqs. (fTTH) and (fTT|) in the BCS and unitarity limits respectively. On the BEG 
side the pairing becomes negligible as compared to the molecular binding energy and 
E ~ \Ei)\/2. The quasiparticle energies are much smaller than the pairing field in the 
center of the trap and are therefore in-gap surface excitations. 

The Bogoliubov-deGennes equations and semiclassical solutions can be extended 
to deformed traps. The quasiparticle energy degeneracy for / = 0, 1, .... will be spht and 
depend on uii. However, due to energy weighted sumrules we expect that the average 
quasiparticle energies can be approximated by Eqs. (|Tn|l and (|TT|l when ojq is replaced by 
This in concordance with the analysis of the resonance superfluidity 



UJ 



[UJlUJ2UJ^ 



1 1/3 



theory which is independent of the symmetry of the cloud ^3] ■ 



4. Radio Frequency Excitation Spectrum 

For detailed calculations of RF response functions for trapped Fermi atoms and at finite 
temperature we refer to recent calculations based on resonance superfluidity theory 
[T3] . Their Bose-Fermion model divides the system of strongly correlated atoms into 
molecular bosons and atoms and introduces effective couplings as well as cutoffs between 
these constituents. They can qualitatively explain the RF response functions at the 
various interaction strengths, densities and temperatures but exceed the effective pairing 
gaps quantitatively. We shall now attempt to understand the RF excitation spectrum 
in terms of the quasiparticle energies calculated above at zero temperature. 

The "effective" pairing gaps reported in |llj are defined at the maximum of the RF 
response function. The excitations energies include the breaking of the pair in states 
|1) + |2) to an in-gap quasiparticle in state |1) and an unpaired quasiparticle in the new 
spin-state |3). The latter energy is, however, subtracted on average as it is related to the 
thermal peak. The RF excitation spectrum is therefore distributed around the lowest 
quasiparticle excitation energies. As discussed above we shall assume that the scale of 
the RF spectral weight in a deformed trap is given by the lowest quasiparticle excitations 
as given Eqs. (jlOj) and (jllll replacing c^o with uj. The distribution is also smeared by 
the spectrum of final h.o. states |3) and by finite temperature effects and damping. The 
spectral response function is a sum over quasiparticle excitation energies weighted with a 
transition matrix elements. These typically decrease rapidly with transition energy [201 
and therefore the maximum of the response function is expected to peak at or near the 
lowest quasiparticle energies. In the BEC limit the maximum of the response function 
is larger than the lowest quasiparticle energy \Ei\/2 by a factor (4/3) only [TTj . 

In Fig. 1 we compare the two lowest quasiparticle excitations Eq and Ei to the 
effective pairing gaps /iz//2 of Ref. [TT] for values of akp around the unitarity limit. 
En are calculated from Eq. (jHl), kp from the central trap densities, and a is calculated 
assuming that the Feshbach resonances resides aX B = 8456*. We observe that the 
effective pairing gaps lie between lowest quasiparticle energies Eq and Ei as argued 
above - except in the dilute BCS limit. The latter departure may be a finite temperature 
effect. In the BCS limit the pairing gap and thus Tc decrease toward the temperature of 
the gas leading to quenching of the gap. Unfortunately, the temperature can presently 
not be measured accurately in traps. 

Theoretically the quasiparticle energies and pairing gaps in units of the Fermi 
energy are universal function of akp around the unitarity limit when np is unchanged. 
The effective pairing gaps at three different Fermi energies obey this scaling to a good 
approximation. It should be noted that the exact position of the Feshbach resonance is 
important otherwise the scaling is destroyed. By placing the resonance in ^Li around 
B ~ 845(5 the effective pairing gaps scale almost perfectly as seen in Fig. 1. Best 
agreement between experimental and calculated collective modes, that will be discussed 
in the following section, is found for a Feshbach resonance residing around B c^ 837G 
For an even lower Feshbach resonance systematic deviations appear between the 
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effective pairing gaps and Ep. 

We emphasize that the pairing gaps are calculated from the semiclassical 
quantization condition of Eq. (jH)) which is based on the Bogoliubov-deGennes equations 
derived from the gap equation in the crossover mode. However, although the crossover 
model provides a simple and qualitatively description of the ground state around the 
BCS-BCS crossover, it differs from the ground state found by Monte Carlo calculations 
and 4-body calculations for both the pairing field and the chemical potential as 
mentioned above. The effect of mean fields were included by replacing the chemical 
potential by that calculated by Monte Carlo. Induced interactions were included by 
using the Gorkov result of Eq. (^ for the pairing field which also is a fair approximation 
around the unitarity limit according to Monte-Carlo calculations |14j . 

5. Collective Modes 

The hydrodynamic and superfluid collective mode frequencies are unfortunately identical 
and we have to rely on other observables to be able to distinguish them and prove, e.g. 
superfluidity. We shall first discuss the collective modes for a very large number of 
particles in a trap and/or with strong interactions. Subsequently, we discuss the effect 
of a finite number of particles in a trap and damping with and without superfiuidity. 

The hydrodynamic and superfiuid collective modes can be calculated in general for 
a polytropic equation of state: P oc n^~^^. Here, the polytropic power is 7 = 1 in a dilute 
interaction dominated BEC, whereas an ideal Bose gas in the normal state has 7 = 2/3 
under adiabatic conditions. A dilute gas of Fermi atoms also has 7 = 2/3 in both the 
hydrodynamic and superfluid limits. Both a Fermi gas ^21 and a BEC j22] has 7 = 2/3 
in the strongly interacting (unitarity) limit. The effective power 7 has been calculated 
at the BCS-BEC crossover for the Leggett model |31 EH] and by Monte Carlo jH] and 
varies between 7 ~ 0.5 — 1.3. 

A spherical symmetric h.o. trap with a polytropic equation of state has collective 
mode frequencies j21] cu^ = UQ{l+2n['j{n+l+l/2) + l\), where / is the angular momentum 
and n the number of radial nodes. In comparison the collective modes in the coUisionless 
limit are those of a free particle: uj/ujq = 2n + I, when its mean free path exceeds the 
size of the cloud. 

In an axial symmetric trap: cji = ti;2 = 00 ± and u^ = Xu!±, the resulting breathing 
modes are the coupled monopole and quadrupole 171 = modes j^Hl- For the very cigar 
shaped traps A ^ 1 used in P and jHH the coupled modes become the radial 



u;rad = V2(7 + l)a;^, (12) 

and axial 



Cc'ax = V3 - (7 + 1)-'CU3 , (13) 

modes. Taking the effective power 7 from the Leggett model at the BCS-BEC crossover 
good agreement is found [211 I2H] with the experiments of jH] and ^U] for the axial and 
and for the radial modes of [H]. The radial mode in ^U] differs and a transition or 
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"break" is observed around the break point B ~ 910G which is accompanied by strong 
damping. 

In the unitarity hmit x = l/lakp) = scahng predicts that 7 = 2/3 as is also 
found in the axial and radial mode of j9.' and in the axial mode of |10j. Furthermore, 
it follows from Eqs. ()12|1 and p3|l that as function oi x = l/{akp) their slopes are: 
^'rad/^i- = 7 v 3/10 and uj'^^/uJ'j, = 7'(3/5)^/^/4 in the unitarity limit. The slope of 7 
and /? can at x = be related as 



The Monte Carlo calculations ^lE], the crossover model and the LOCV model [T^ IT3] 
all find f3'{x = 0) ~ —0.10 and (3{x = 0) ~ —0.56 (the crossover model does, however, 
give I3{x = 0) = —0.42 which probably is due to the lack of Hartree-Fock energy 
corrections). Thus we obtain 7'(x = 0) ~ 0.40. The resulting slopes of the collective 
frequencies u'^^^ and tu^^ are also compatible with experiment [HI HD] in the unitarity 
limit. 

The measured damping of the modes is not compatible with hydrodynamics. As 
pointed out in jH] the damping rate increase with increasing temperature whereas the 
opposite is expected in hydrodynamics. In a superfiuid the condensate is gradually 
depleted as the temperature increases and coupling between the normal and superfiuid 
components increase damping as observed. The damping in ^Uj peaks at the transition 
which together with the abrupt transition in frequency indicates a superfiuid to 
coUisionless transition rather than a smooth hydrodynamic to coUisionless transition. 

The damping rate from coUisional damping in a normal gas can be estimated from 
the semi-quantitative form for the transition between coUisionless and hydrodynamic 
hmit 123 EOl 

u' = ul- ^^-^^ . (15) 

1 - lUJT^oll 



For the radial mode the hydrodynamic frequency is ujh — a/IO/Scu^ whereas the 
coUisionless is uc — 2uj^. The relaxation processes are expressed in terms of a collision 
time Tcou- The maximal damping of the collective mode Im{uj)/uj occurs between the 
hydrodynamic and coUisionless limits for ujTcou = 1 and is: Im{uj) ^ 0.09co'x, for the 
radial mode. The damping rates in ^Oj and J2Hj are considerably larger around the 
break points, and can thus not be caused by coUisional damping alone. 

Alternatively, the enhanced damping and break point could be caused by a 
superfiuid to coUisionless transition. However, in the center of the trap the critical 
temperature T^ = 0.28Epexp{TTx/2) is at the break point x ~ —0.5 much larger than 
the gas temperatures, T;^0.03Ef- Thus only the low density surface layer is not in 
the superfiuid phase. Estimates of the effective scattering cross sections indicate that 
the gas is coUisionless in the normal phase at the very low temperatures present in the 
experiments of JJ]] and [01 EH] • 

It was pointed out in Ref. jTU] that the lowest quasiparticle excitation energy is 
comparable to the collective energy hLj± of the radial mode at their transition. An 
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interesting coupling between the single particle states and the collective mode may 
therefore take place which is special for a finite system. This coupling can be studied in 
the RPA equations which describe how the collective modes are build up of quasiparticle 
states. The collective modes are calculated for dilute and spherical h.o. traps in J2ZI by 
solving the RPA equations using the quasiparticle states and energies as input. When 
interactions are weak such that 2Gln(e'"ni?)^/ia;o, the collective modes are dominated 
by the pairing gap and the h.o. shell structure. For stronger interactions the collective 
modes can be calculated from RPA and they approach those of a superfluid. 

Though the collective modes for both the dilute system and the unitarity regime 
typically are of order the harmonic oscillator frequency, the underlying quasiparticle 
spectrum is very different. In the dilute limit the pairing gap and the lowest excitation 
energies can be far below collective energies. This is similar to the situation in atomic 
nuclei where pairing energies are of order IMeV whereas giant resonances lie around 
10-20MeV. In contrast, for trapped atoms in the unitarity regime the quasiparticle 
excitation energies En lie well above hu. Towards the BCS limit, however, Eq approaches 
hurad- It follows froui the RPA equations that the collective frequency both decreases 
and is damped with respect to a system with an infinite number of particles. This is in 
qualitative agreement with the experiment of Ref. jTUj on the BCS side. 

The observed transition or break in the radial frequency occurs at B = 910G ^UI for 
both i'± = 750Hz and 2.4kHz, i.e. for the same scattering length but for two different 
kp and radial frequencies. However, we find that in both cases the resulting radial 
frequency and lowest quasiparticle energy almost coincide at the break point, i.e. 

hUrad ^ 2Eo . (16) 

Note that Ep and kp do not scale linearly with i>± as does Wrad- It is thus a coincidence 
that Eq. ()l(i|l is fulfilled for the same scattering length and therefore at the same B 
for the two different i'±. In the experiment of Ref. [H] both the trap frequencies are 
larger than in ^Hj and one can calculate that the break should occur around B ~ 980G 
according to the condition in Eq. (fTI)|) . Preliminary data J2H] finds that damping does 
increase at such magnetic fields but that a possible break point is situated at even higher 
magnetic fields. 

More precise measurements of Tc and T are required before the precise phase, 
damping mechanisms and break points can be determined. 

6. Summary and Outlook 

We have argued that the recent experiments on RF spectroscopy |.n| and collective 
modes [HI EIj give strong evidence for superfluidity in traps with Fermi atoms. The 
results compare very well with theoretical calculations [13 1211 121] and the quasiparticle 
excitation energies discussed above as function of interaction strength, density and 
temperature. 

There are, however, some details and deviations between theory and experiments 
that needs further investigation. The effective pairing gaps of Ref. ^I] only partially 
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obey the scaling with Ep. The radial mode collective frequency of Ref. ^U] violate 
the scaling result in the unitarity limit and undergo a transition. Better control of the 
temperature and the exact position of the Feshbach resonance as well as experiments 
at several densities and temperatures are all desirable. These are necessary for a 
detailed check of the scaling with akp around the unitarity limit and to pin down 
transitions and critical temperatures between superfluid and normal phases, coUisionless 
and hydrodynamic dynamics. 

In this respect we can appreciate the great advantage of atomic gases, namely 
the great number of tunable parameters such as the number of particles, densities, 
interaction strengths, temperatures, trap deformation, number of spin states, etc. They 
therefore hold great promise for a more general understanding of pairing phenomena 
atomic Fermi gases but also in solids, metallic clusters, grains, nuclei, neutron stars, 
quark matter, etc. 

Acknowledgements: Discussions with J. Thomas and R. Grimm are gratefully 
acknowledged. 
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